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Abstract 

We study constraints on time variation of the fine structure constant a from cosmic microwave 
background (CMB) taking into account simultaneous change in a and the electron mass nie which 
might be implied in unification theories. We obtain the constraints —0.097 < Aa/a < 0.034 at 
95% C.L. using WMAP data only, and —0.042 < Aa/a < 0.026 combining with the constraint on 
the Hubble parameter by the HST Hubble Key Project. These are improved by 15% compared 
with constraints assuming only a varies. We discuss other relations between variations in a and 
nie but we do not find evidence for varying a. 
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I. INTRODUCTION 



One of fundamental questions in physics is whether or not the physical constants are 
literally constant. In fact, the physical constants may change in spacetime within the context 
of some unification theories such as superstring theory and investigating their constancy 
is an important probe of those theories. In addition to the theoretical possibility, some 
observations suggest the time variation of the fine structure constant (the coupling constant 
of electromagnetic interaction) a. Recent constraints including such non-null results are 
briefly summarized as follows. 



Terrestrial limits on a come from atomic clocks 
reactor in Gabon 



the Oklo natural fission 



and meteorites. Ref. derived the limit on the temporal derivative of a 
at present as (—0.3 ±2.0) x lO^^^yr^^. Measurements of Sm isotopes in the Oklo provide two 
rounds on the variation of a as Aa/a = — (0.8±1.0) x 10~^ and Aa/a = (0.88±0.07) x 10^''' 
^. The former result is null, however, the latter is a strong detection that a was larger at 
z ~ 0.1. The meteorite bound obtained by measuririg ^^^Re decay rate is now controversial 
although varying a is not suggested anyway. Ref. jy] obtained Aa/a = (—8 ± 8) x 10"'', 



whereas Refs. 



10 



111 ] argued that the constraint should be much weaker due to uncertainties 



in the decay rate modeling. 

On the other hand, there are three kinds of celestial probes. One is to use big bang 



nucleosynthesis (BBN) 



12 



13 



15l | which provides constraints at very high redshifts 



10^-10^°), for example, -5.0 x lO'^ < A«/« < 1.0 x 10"^ (95% C.L.) (jj or |Aa/al< 



6 X 10-2 



19 



i4. 



21 



22 



'he second is from the spectra of high-redshift quasars {z ~ 1-3) |lj 



where there are conflicting results. Refs. [17|, L 
was smaller at z > 1, l^aja = (—0.543 ± 0.116) x 10"^ Q]. This result, however, is 

. [19| obtained 



2C 



21 



22 



suggested that a 
3sul 

. For example, Ref. 



not supported by other observations 3, 
the constraint as Aa/a = (—0.6 ± 0.6) x 10"^ and the others too found no evidence for 

varying Finallv, we can use Cosmic Microwave Background (CMB) to measure a at 

Q, 3, Q- Analyses using pre-WMAP data are found in Refs. Q, 22, Q- 



z ~ 1100 12 

3. 



Refs. 



30| derive a constraint using the WMAP first-year data, —0.05 < Aa/a < 0.02 



or —0.06 < Aa/a < 0.01 (95% C.L.) respectively with or without marginalizing over the 
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running of the spectral index js^. ^ 

Recently, there are many studies on constraining the time variation in a which accompa- 
nies the variation of the other coupling constants, as would occur rather naturall y in unified 



theories. Under such a framework, BBN has been studied in Refs 



40 



in Ref. 



meteorites in Ref. 



43, 



Ic 35 



37|, 



38, 



3S 



the Oklo reactor 



4J|. For example, BBN constraint 



^ quasar absorption systems in Ref. 
I, and atomic clock experiments in Refs. 
improves by up to about 2 orders (—6.0 x 10""^ < Aa/a < 1.5 x 10~^ |l3]) although the fac- 
tor may vary depending on how they are correlated to a. Hence, it is important to consider 
the variation of other coupling constants along with a. 

In this paper, we investigate constraint on the time variation of a from CMB using the 
WMAP first-year data. In particular, we consider mg to vary dependently to the variation 
in a because, as mentioned above, such might be the case in some unified theory. Since such 
theory is now under development, how their variations are related to each other can not be 
predicted. Therefore, we work in a phenomenological way guided by a low energy effective 
theory of a string theory and adopt to vary vrif. in power law of a. ^ 

In the next section, we briefiy review the recombination process in the early universe and 
make clear how it depends on a and nif.. In section IIIIl we illustrate their effects on the 
epoch of recombination and the shape of CMB power spectrum. In section HVl we describe 
the relations between the variations in a and mg which we adopt in our analysis. In section 
IVl we present our constraint and we conclude in section IVTl 



II. a AND me DEPENDENCE OF THE RECOMBINATION PROCESS 

Non-standard values of a and me modify the CMB angular power spectrum mainly by 
changing the epoch of recombination. Thus, let us visit briefly the recombination process 
in the universe and see where those constants appear. We follow the treatment of Ref. ^^k, , 



^ Note that two results suggesting varying a, one of the Oklo results y and the quasar observation of 
Ref. Ql have different signs of Aa. These results, if correct, can not be explained by homogeneous and 
monotonically time- varying a. They may indicate that a is not a monotonically varying function of time 
or as investigated in Refs. |3l[l33.l3^ . may suggest a spatial variation of a. In passing, we refer to Ref. [s^ 
for th p pfl Fer. t s of spatial varying a on CMB. 

^ Refs. |45|,|4g studied CMB constraints on the variation of Higgs expectation value whose effect is assumed 
to appear only in the change in rUe so there is some overlap between their analysis and ours. However, 
they did not consider a variation at the same time. 



3 



which is implemented in the RECFAST code. They have shown that the recombination 
process is well approximated by the evolutions of three variables: the proton fraction Xp, 
the singly ionized helium fraction xhcii, and the matter temperature Tm- Their equations 
are given bellow. We denote the Boltzmann constant k, the Planck constant h and the 
speed of light c. In addition to the variables above, Xp = Up/nu and xhcii = nueii/nu, we 
use the electron fraction Xe = Ue/nn = Xp + xhcII as an auxiliary variable (note that nx 
stands for the number density of species X but nn is defined as the total hydrogen number 
density, including both protons and hydrogen atoms), z is used for the redshift and H for 
the expansion rate. 

Adopting the three level approximation, the time evolution of the proton fraction of Xp 
is described by 



'2^ _ 

dz 



[xeXpUnRn - ^h(1 - Xp)e-'''^«/'=^*^] , 



H{z){l + z) 

(1) 

where i^h = c/(121.5682 nm) is the Lya frequency and i?H is the case B recombination 
coefficient for H which is well fitted by 

i?H = lO-i^F^^ m^s-i (2) 



with t = Tm/(10^ K), a = 4.309, b = -0.6166, c = 0.6703, d = 0.5300 j4S|, and the fudge 
factor F = 1.14 introduced to reproduce the more precise multi-level calculation [4^. is 
the photoionization coefficient 

3 

and Ch is the so-called Peebles reduction factor 

^ ^ [1 + KuAnriujl - Xp)] 
" [1 + Kn{Au + f3u)nn{l - Xp)y ^' 

where the binding energy in the 2s energy level is -Bh2s = 3.4 eV, the two-photon decay rate 
is Ah = 8.22458 s'^ and Kh = c^/{Stiu^H). 

Here, comments on what Eq. (Q) means may be in order. The first term in the square 
brackets in Eq. (Q) represents the recombinations to excited states of the atom, ignoring 



recombination direct to the ground state. The second term represents the rate of ionization 
from excited states of the atom. The difference of those two terms is the net rate of produc- 
tion of hydrogen atoms when one could ignore the Lya resonance photons. These photons 
reduce the rate by the factor Ch, which can be written as the ratio of the net decay rate to 
the sum of the decay and ionization rates from the n = 2 level, 

Ar + Ah 
" Ak + Ah + /5h- 

We have rewritten Eq. (jlj) to derive this expression using the decay rate Ar = (i^n'^is)"^ 
allowed by redshifting of Lya photons out of the line, and using riis = nn — Up which is 
justified by the far greater occupation number of the hydrogen atom ground state than that 
of the excited states altogether. 

The evolution of the singly ionized helium fraction of xhgii is similarly described by 

dz ~ H{z){l + z) 

(6) 

where /hg is the total number fraction of helium to hydrogen (using primordial helium mass 
fraction Yp, fue = Yp/{4:{1 - Yp)} where we take Yp to be 0.24), z/Rei = c/(60.1404 nm) is 
the frequency corresponding to energy between ground state and 2^s state, and i?Hei is the 
case B Hel recombination coefficient for singlets js^ 



I-P / rp \ 1+P ^ 



m^s \ 



(7) 



with g = 10 16 '^^^, p = 0.711, Ti = 10^'^^^K and T2 fixed arbitrary 3K. jSuei is the photoion- 
ization coefficient 



/3HeI = i?HeI J exp ^-^J , (8) 

Cneii = 

[1 + /^HelAHenH(/He - a^Hell) exp( AE/fcTAf )] 
[1 + i^Hel(AHe + /?Hel)^H(/He " a;HeIl) exp{AE/kTM)] ' 

(9) 



where the binding energy in the 2s energy level is -BneKs = 3.97 eV and the two-photon decay 
rate is Ahc = 51.3 s~^, -f^Hei = c^/(8^^Hei-^) the energy separation between 2^s 

and 2^p, AE/h = c/ (58.4334 nm) — z/Hei. Note that, contrary to the case with H, the energy 
separation AE between 2^s and 2^p is so large that we can not neglect it [sO^ . 
The matter temperature Tm is evolved as 

dTM 



dz 



SaraRTn Xe „ . 2Tm 

■[-LM — -LR) + 



3H{z){l + z)m,l + Me + Xe' (1 + ^) 



(10) 



where Tr is the radiation temperature, = 2a'^h'^/{3'Kmlc'^) is the Thomson cross section, 
and = /{1207Tc^h^) is the black-body constant. The Compton scattering makes Tr and 
Tm identical at high redshifts. However, the adiabatic cooling becomes dominant at low 
redshifts, which leads to the significant difference between Tm and Tr. 

Now, we explain how quantities which appear in these equations depend on a and rUe. 



Two-photon decay rates scale as a^rrie 



53|. Since binding energies scale as a'^rrie, so do 



i^H and z/Hei, and and K^^i scale as a~^m~^. The remaining task is to investigate how 
the recombination coefficient R depends on a and mg. To do this, we follow the treatment 
of Ref. [2^. The recombination coefficient can be expressed as 



«^i:i-(-^^)i(|^)'-(|-) 



n,l 

X r ^f!^, (11) 

JB^/kTM ^My) - 1 

where -B„ is the binding energy for n-th excited state and (T„; is the ionization cross section 
for (n,Z) excited state. The asterisk in the upper bound of summation indicates that the 
sum needs to be re gula ted, but since this regularization depends only weakly on a and rrie, 
it can be neglected |0] . The cross section ani scales as a~^m~'^ [53]. Altogether, 

dR{TM) 2 / _ , ^ dR{TM 



R{Tm)-Tm^^], (12) 
oa a \ oTm J 

dRiTu) _ 1 / dR{TM) \ 

- K^^^^'^ + ^^^f^ ) • ^^^^ 

Combining with the fitting formulae (0) and ((Tj), we obtain how i?H and -Rnei depend on a 



and m, 



e- 




FIG. 1: The ionization fraction Xe as a function of redshift z for no change of a (sohd curve), an 
increase of a by 5% (dashed curve), a decrease of a by 5% (dotted curve). 




FIG. 2: The ionization fraction Xg as a function of redshift z for no change of rrie (sohd curve), an 
increase of nie by 5% (dashed curve), a decrease of rUe by 5% (dotted curve). 

III. EFFECTS ON THE EPOCH OF RECOMBINATION AND CMB ANGULAR 
POWER SPECTRUM 



We have investigated the equations which describe the process of recombination and 
how they depend on the couphng constants in the previous section. We incorporate the 
dependence on a and rrie into the RECFAST code and solve the equations for the 
ionization fraction Xe as a function of redshift with several different values of a and rrie- The 
results are shown in Figs. Q and El We have assumed a flat universe and used cosmological 
parameters {ub,Um,h) = (0.024,0.14,0.72), where Ub = flbfi^ is the baryon density, = 
VLmh? is the matter density, h is the Hubble parameter, and VL denotes the energy density 
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FIG. 3: The visibility function as a function of conformal time for no change of a (solid curve), an 
increase of a by 5% (dashed curve), a decrease of a by 5% (dotted curve). 




Ame/me- +0.05 
Ame/me= 
Arrfe/me=-0.05 
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FIG. 4: The visibility function as a function of conformal time for no change of me (solid curve), 
an increase of mg by 5% (dashed curve), a decrease of me by 5% (dotted curve). 

in unit of the critical density. The most important feature is the shift of the epoch of 
recombination to higher z as a or me increases. We can also see this by the rightward shift 
of the peak of the visibility function shown in Figs. El and IH This is easy to understand 
because the binding energy scales as o^vUf. and photons should have higher energy to 
ionize hydrogens. 

In Figs. El and El we show the power spectrum of the CMB temperature anisotropy for 
several different values of a and me as calculated by the CMBFAST code with the 
modified RECFAST. We consider a flat ACDM universe with a power-law adiabatic pri- 
mordial fluctuation. The adopted cosmological parameter values are {uJb,u)m, h,ns,T) = 
(0.024,0.14,0.72,0.99,0.166) where r is the reionization optical depth and is the scalar 
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V 

FIG. 5: The spectrum of CMB fluctuations for no change of a (sohd curve), an increase of a by 
5% (dashed curve), a decrease of a by 5% (dotted curve). 
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FIG. 6: The spectrum of CMB fluctuations for no change of nie (solid curve), an increase of me 
by 5% (dashed curve), a decrease of mg by 5% (dotted curve). 

spectral index. We fix the value of the amplitude of primordial power spectrum. We can see 
two effects of varying a and mg in Figs. El and IHl Increasing a or rrie shift the peak positions 
to higher values of I and amplify the peak heights. 

The peak position shift is understood as follows. Using Ip to denote the position of a 
peak, r0{z) for the angular diameter distance and rs{z) for the sound horizon, one can write 



Ip ~ 



where zis is the redshift of the last scattering surface. Increasing a or rrie increases the 
redshift of the last scattering surface due to the larger binding energy, as in Figs. El and El 
The higher zis in turn corresponds to a smaller sound horizon and a larger angular diameter 
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distance, which lead to a higher value of Ip. 

The changes in the peak heights are caused by modifications to the early ISW effect 
and the diffusion damping. The larger zis leads to the larger early ISW effect making the 
first peak higher. To consider the effect beyond the first peak, we focus our attention on 
the visibility function. The peak of the visibility function moves to a larger redshift since 
the recombination occurs at higher redshift, when the expansion rate is faster. Hence, the 
temperature and Xe decreases more rapidly, making the peak width of the visibility function 
narrower (see Figs. El and 0]). Since the width of the visibility function corresponds to the 
damping scale, an increase in a or rrie decreases the effect of damping. This is the reason 
why the amplitude at larger / increases with increasing a and mg. Moreover, as seen in 
Figs. inland 01 a changes the visibility function width more than rrie does (quantitatively, an 
increase of a or rrie by 5% makes the full width at half maximum of the visibility function 
narrower by 10% or 2% respectively) because the binding energy which scales as a^mg. Thus 
the damping scale is more sensitive to the change of a than rrie, as appears in Figs. El and El 

Now we discuss the effects of varying a and rrie somewhat more quantitativelv using 
the following four quantities which characterize a shape of CMB power spectrum j^: the 
position of the first peak li, the height of the first peak relative to the large angular-scale 
amplitude evaluated at / = 10, 

H.^i^y (15) 



vATio 

the ratio of the second peak (/2) height to the first 



the ratio of the third peak (/s) height to the first 



AT,. 



^^.-1^1 (16) 



where (AT^)^ = /(/ + l)C//27r. Note that these four quantities do not depend on overall 
amplitude. We calculate the response of these four quantities when we vary the parameters 
cjb, ^ni) h, T, Us, a and rrie- When we vary one parameter, the other parameters are 
fixed and fiatness is always assumed (especially, increasing h means increasing because 
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Ah = 16^ - 25^ - 47^ 

An, Aa Amp , , 

+36 ^ + 290 + 150 (18) 

Us OL rrie 

AH, = 3.0^ - 3.0^ - 2.2^ - 1.7^ 

Aus Aa , Ame , 
+18 ^ + 3.9 + 1.4 (19 

Us OL me 

Aujh Aujm An, 
AHo = -0.30 ^ + 0.015 ^ + 0.41- ' 



LOb Um Us 

+0.91^ + 0.30^, (20) 

a me 

AH, = -0.19^ + 0.21^ + 0.56-^^' 



tOb ris 

+0.57^-0.019^, (21) 
a me 

and values at the fiducial parameter values are h = 220, H, = 6.65, H2 = 0.442 and 
H, = 0.449. Derivatives of li and Hi with respect to a and m^ are positive as is expected 
from the considerations above. Furthermore, A/i/Aa and A/i/Amg are much larger than 
the other derivatives of li while AHi/Aa and AH,/ Ame have relatively similar values to 
the other derivatives of H,. Since such changes are most effectively mimicked by the change 
in /i, it is considered to be the most degenerate parameter with a and me- We have seen 
above that when a or me increases, the first peak is enhanced by larger ISW effect and the 
second or higher peaks are enhanced by smaller diffusion damping. The derivatives of H2 and 
H, tell us which effect is important. Since AH2I Aa and AH,/Aa are positive and larger 
than the derivatives with respect to mg, we see that the effect on the diffusion damping 
is more significant than that on the early ISW for varying a. They seem to somewhat 
cancel each other for varying me especially regarding H,. Such behavior is consistent with 
the consideration at the end of the previous paragraph, that the diffusion damping is more 
sensitive to the change in a than me- 

IV. RELATION BETWEEN VARIATIONS OF a AND rUe 

We expect a unified theory can predict the values of the couphng constants, how they are 
related to each other and how much they vary in cosmological time scale. In string theory, a 



11 



candidate for unified theory, there is a dilaton field whose expectation value determines the 
values of coupling constants. However, since it is not fully formulated at present, we have 
to assume how a and vHf. are related to vary and constrain their variations. To be concrete, 
following Ref. js^, let us start from considering the low energy action derived from heterotic 
string theory in the Einstein frame. The action is written as 



S = I d'^x^ ( —R - 



a' 



- ij^^D''^ - n-'m^^^ - j^n^F^^F'^'^j , (22) 

where $ is the dilaton field, (p is an arbitrary scalar field, and is an arbitrary fermion. 
Dfj^ is the gauge covariant derivative corresponding to gauge fields with field strength F^,^, 
= SvrG and Q = e"*^*/^ is the conformal factor which is used to move from string frame. 
More concretely, is the Higgs field and V{(f)) is its potential. The overall factor Q before 
the scalar potential means that the Higgs vacuum expectation value (H) is independent of 
the dilaton so it is taken to be constant. Ff^i, is the gauge field with gauge group including 
SU(3)xSU(2) xU(l). We define its Lagrangian density for the gauge field as —{l/Ag'^)F^yF^^'' 
where g is the unified coupling constant. Compared with equation Eq. (j221), 

1 a'e"^"* 
'^{I^ ^ Ak^ ' ^^^^ 

where Mp is the Planck scale. We can calculate the gauge coupling constants at low energy 
using renormalization group equations, a almost does not run, and hence the a at low 
energy 

a ^ a = = —. 24 

47r TTa 

As for the other gauge coupling constants, variation of the strong coupling constant may 
affect CMB since its low energy value determines the QCD scale Aqcd which in turn deter- 
mines nucleon masses. However, how the variation of Aqcd is related to that of a at low 



energy can not unique 
of unification scheme 



y be determined from eq. and especially depends on the details 



57[. Therefore, for simplicity, we just assume Aqcd does not vary. 



The ip^s are the ordinary standard model leptons and quarks. As we take (H) = const, the 
Yukawa couplings depend on the dilaton as e"*/^. Therefore the relation between variations 



12 



of a and mg is given by 

1 /2 

me + Ame _ / 01 + Aa \ 
me V « / 

In this paper, we also consider other possibihties phenomenologically by adopting a power 
law relation as 

me + Ame _ + Aa Y 
me V « / ' 

and compute constraints for several values of p. In addition to the case with changing only 

a = 0) and the model described above (p = 1/2), we consider cases with p = 2 and 4. 



V. CONSTRAINTS ON VARYING a AND me 



We constrain the variation of a in the models described in the previous section using 
the WMAP first-year data. CMB power spectrum is calculated by CMBFAST jssj with 



RECFAST 



471 ] modified as in Sec. El The is coraputed for TT and TE data set by the 



likelihood code supplied by the WMAP team 



58, 



59 



6011. As defined in Sec. IIIIl we consider 



six cosmological parameters Ub, ujm, h, r, Ug and overall amplitude A in the ACDM model 
assuming the flatness of the universe. We report A in terms of /(/ + l)Ci/27r at / = 2 in unit 
of fiK^. In this paper we do not consider gravity waves, running of the spectral index and 
isocurvature modes. We calculate minimum as a function of a and derive constraints on 
a. The minimization over six other parameters are performed by iterative applications of 



the Brent method 



of the successive parabolic interpolation. More detailed description 



of this minimization method is found in Ref. ^ . We search for minimum in the region 



T < 0.3, which is a prior adopted in Refs. 2^ 



We derive constraints with or without 



the constraint on the hubble parameter h. When we combine it, we use the Hubble Space 



Telescope (HST) Hubble Key Project value h = 0.72 ± 0.08 |63[ whose error is regarded as 
gaussian la. 

Fig. 13 shows our results of minimization. It compares varying a only and varying a 
and me with the relation of eq. ()25|1 . respectively with or without the HST prior. Without 
the HST prior, we find at 95% CL. that —0.107 < Aa/a < 0.043 with changing a only 
and —0.097 < Aa/a < 0.034 with the model described in the previous section. Although 
the best fit a is 4% less than the present value, we find that Aa = is consistent with the 
WMAP observation and evidence for varying a is not obtained. The effect of varying me 
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- a only a only with HST prior 

- a and me a and nig witii HST prior 




FIG. 7: function for variation in a only (solid curves) and a and mg simultaneously with the 
relation of eq. (dashed curves) . The HST prior is imposed for the thinner curves. 




FIG. 8: function for variation in a and rrie with several values of p defined in Eq. (|26|) . The 
curves are for, from outside to inside, p = 0, 1/2, 2 and 4. The HST prior is imposed. 

simultaneously is found to make the constraint more stringent by 13%. This rather small 
effect is reasonable since, as is discussed in Sec. lllH the effect of me on CMB power spectrum 
is slightly smaller than a, and the relation of eq. (j25|) we adopt here does not change me 
much relative to a. 

We find that minimum is given at Aa/a = —0.04 with {ujb,ujm, h,ns,T, A) = 
(0.021,0.132,0.523,0.979,0.146,942) for the case of changing a only, and Aa/a = -0.04 
with {uJb,uJm,h,ns,T,A) = (0.020,0.131,0.485,0.979,0.140,907) for the case of changing a 
and me together. Both cases have notably small values of h. Since h is considered to be the 
most degenerate parameter with a or me as discussed in the end of Sec. IllH it is instructive 
to investigate how constraints tighten when h is limited to higher values such as the HST 
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power law index 


constraint (95% C.L.) 


p = 


-0.048 < Aa/a < 0.032 


p = l/2 


-0.042 < Aa/a < 0.026 


p = 2 


-0.031 < Aa/a < 0.017 


p = 4 


-0.023 < Aa/a < 0.011 



TABLE I: The constraints on a with mg varies according to various power law relation eq. (|26l) of 
the index p. The HST prior is imposed. 

measurement. From Fig. [71 we obtain, with the HST prior, that —0.048 < Aa/a < 0.032 
with changing a only, and —0.042 < Aa/a < 0.026 with the model described in the previous 
section. Compared with no HST prior constraints, they are stringent by about factor of 2 
for both cases. Moreover, since low values of h which give good fit with Aa/a ^ —0.04 are 
ruled out by the HST prior, the center of allowed region has shifted to larger Aa. 

Here, we comment on ..e con.... p.ev,o.„ oUa.ed ,y Re.. HQ t.e WMAP 
data. As mentioned in SeclU they reported the constraint on a to be —0.06 < Aa/a < 0.01 
(95% C.L.). They fixed rrie when varying a and values quoted here is the case with no 
running for the primordial power spectrum. This constraint seems to have been obtained 
with marginalization on grid with < Q\ < 0.95 2^ so it should be compared with our 
constraint without the HST prior, —0.107 < Aa/a < 0.043, which is much weaker than 
theirs. The difference might be traced to the different analysis method but we could not 
reproduce their results by our method. 

Finally, we investigate the cases in which rrie varies more than a. We consider the models 
with p = 2 and 4 in eq. (f^H|) . We calculate constraints with the HST prior and results are 
summarized in Fig. |Hland Table H] along with p = and 1/2 cases. Compared with p = 
(only varying a) case, the constraints become smaller by 40% {p = 2) and 60% (p = 4). 
Although those constraints are much smaller than the case with p = 1/2, they are still 
consistent with Aa = 0. 
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VI. CONCLUSION 



In summary, we have studied a CMB constraint on the time varying fine structure con- 
stant a taking into account simultaneous change of electron mass me which might be implied 
in superstring theories. We have searched sufficiently wide ranges of the cosmological pa- 
rameters and obtained the WMAP only constraint at 95% CL. as —0.097 < Aa/a < 0.034. 
Combining with the measurement of h by the HST Hubble Key Project, we have obtained 
more stringent constraint as —0.042 < Aa/a < 0.026, which improvement is explained by 
the strong degeneracy between a or rrie and h. These constraints, obtained by adopting 
the model with eq. (|25|) are only slightly tighter than those assuming only a variation: 
-0.107 < Aa/a < 0.043 and -0.048 < Aa/a < 0.032, without and with the HST prior 
respectively. This is reasonable since the effect of rrie on CMB power spectrum is similar to 
that of a (sec. IIII|) and we adopted to vary rrie milder than a ( eq. (j2Sl))- 

We have also considered other possibilities for the relation between Aa and Amg as in the 
form of eq. (j26|l with p = 2 and 4. In these cases, the constraints become tighter by roughly 
a factor of two. This may not look as drastic as the BBN bounds but CMB bounds are 
promising and have advantage that there will be future experiments with higher sensitivities 



to a as investigated in Refs. 29|, |30|. In this paper, we do not find evidence of varying a 



in the CMB data of WMAP. We will wait and see whether future experiments give us more 
stringent bound on a or evidence for varying a. 
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